OCU-PHYS-205 



Perturbative out of equilibrium quantum field theory 
beyond the gradient approximation and generalized 

Boltzmann equation. 

H. Ozaki* 

Graduate school of science, Osaka City University, 
Sumiyoshi-ku, Osaka 558-8585, Japan 

(November, 2003) 



Abstract 

Using the closed-time-path formahsm, we construct perturbative frameworks, in 
terms of quasiparticle picture, for studying quasiuniform relativistic quantum field 
systems near equilibrium and non-equilibrium quasistationary systems. We employ 
the derivative expansion and take in up to the second-order term, i.e., one-order higher 
than the gradient approximation. After constructing self-energy resummed propagator, 
we formulated two kind of mutually equivalent perturbative frameworks: The first one 
is formulated on the basis of the "bare" number density function, and the second one 
is formulated on the basis of "physical" number density function. In the course of 
construction of the second framework, the generalized Boltzmann equations directly 
come out, which describe the evolution of the system. 



1 Introduction 

Ultrarelativistic heavy-ion collision experiments have begun at the BNL Relativistic Heavy 
Ion Collider (RHIC) and will soon start at the CERN Large-Hadron Collider (LHC) in 
expectation of producing quark-gluon plasma (QGP) A produced system in such a 

experiment is nonequilibrium and evolves into confined (or hadronic) and chiral-symmetry 
breaking phase through phase transitions. For dealing with nonequilibrium quantum-field 
systems, several formalisms are available. Among those are the closed-time-path (CTP) 
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formalism OE] and the nonequilibrium thermo field dynamics In this paper, we employ 
the former. 

Throughout this paper, we are interested in quasiuniform relativistic quantum field sys- 
tems near equilibrium and nonequilibrium quasistationary systems, which we simply refer 
to as "out of equilibrium" systems. The out of equilibrium systems are characterized by two 
different spacetime scales jlj. The first one is microscopic or quantum- field theoretical scale, 
which characterizes the range of radiative corrections to the reactions taking place in the 
system. The second scale is macroscopic or statistical scale, which measures the relaxation 
of the system. For a weak coupling theory, in which we are interested, the former scale is 
much smaller than the latter one. 

The "two-scale structure" mentioned above is implemented as follows. Let G{x, y) be 
a generic two-point function. G{x,y), with x — y fixed, does not change appreciably in 
{x + y)/2. The relative coordinates x — y are chosen as the microscopic coordinates while 
X = {x + y)/2 are chosen as the macroscopic coordinates. A Wigner transformation, i.e., a 
Fourier transform with respect to x — y, yields 



[P^ = p)]. Then in Eq.(^ can be regarded as the momentum of quasiparticle partic- 
ipating in the microscopic reaction under consideration. Because of the weak X-dependence 
of G{X; P), the derivative expansion makes sense. 



where dyt^ = d/dY^. We refer to the first term, G{Y\P), on the right-hand side (RHS) 
as the leading or zeroth-order part (term), to the second term to the gradient or first-order 
part (term), and to the third term as the second-order part (term). Most works reported as 
far ( see, e.g., jEllZlE]) on out-of-equilibrium perturbation theory are carried out within the 
gradient approximation, i.e., the approximation of taking up to the gradient term in Eq.(|^. 

The purpose of this paper is to frame perturbative frameworks beyond the gradient 
approximation by taking up to the second-order term on the RHS of Eq.©. We construct 
two mutually equivalent perturbative frameworks. One framework is constructed in terms 
of the "bare" -number density function, and the other, which is called physical- iV scheme, is 
constructed in terms of the "renormalized" -number density function. For definiteness, we 
take up an out-of-equilibrium complex-scalar quantum field system^. Generalization to other 
theories is important. In this relation, it has been pointed out that a care should be taken [H] 
when one applies the derivative expansion to gauge theories. 

^This system has been dealt with in [7] to the gradient approximation. Fermion and gauge-boson systems 
have been treated in [H]. 
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The perturbative framework to be constructed accompanies the generahzed Bohzmann 
equation (GBE) for the number density of quasiparticles. The framework allows us to com- 
pute any reaction rate by using the reaction rate formula jH]. Substituting the computed 
net production rates of quasiparticles into the GBE, one can determine the number densities 
as functions of macroscopic spacetime coordinates X'^, which describes the evolution of the 
system. 

The plan of the paper is as follows. In Sec. |2l on the basis of the CTP formalism supple- 
mented with a quasiparticle representation, we construct the free quasiparticle propagator 
up to the second-order term in the derivative expansion, Eq.Q. In Sec. El we construct 
the self-energy resummed propagator and observe that large contributions come out in gen- 
eral. In Sec. m after introducing "physical" number density, we briefly describe "the bare-iV 
scheme." In Sec. imposing the condition that no large contributions appears, we construct 
a "healthy" perturbative framework (the physical- scheme). It is shown that, on the en- 
ergy shell, the condition turns out to the generalized Boltzmann equation. Furthermore we 
impose the condition that, on the energy shell (with p = 0), the medium effect as well as 
the quantum correction to the "mass" vanish, which leads to a gap equation. In Sec. IHl 
the similarity between the bare-iV scheme and the physical- iV scheme is discussed. Sec. dis 
devoted to summary. 

2 Closed-time-path Formalism 
2.1 Preliminaries 

The Lagrangian (density) of the complex-scalar field theory reads 

C = -<p\x) {dl + m') 0(x) - ^ (0t(a;)0(x))' . (3) 

To obtain an efficient perturbative scheme, we introduce a weakly x-dependent mass function 
M{x, y) pHI and write C as 

(4) 

)] (5) 

</>(y), (6) 

(7) 
where 

M\x,y) = e{i{d,o-dyo))Ml{x) + e{-i{d,o-dyo))Ml{x) (8) 



c 

Co 
r 



Co + Cci + £int5 

-0t(x) J [5\x - y)dl + M^{x, y 
(j)^{x) j d'^y A4'^{x,y) — m^5{x — y) 

~4 
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d'P g-.P-(.-.)^2^^. (X = (x + y)/2) , (9) 



(2vr) 

A<2(X;po) = ^(po)M2(X) + 0(-Po)M!(X). (10) 

Here M_|_(M_) is the effective mass of a quasiparticle (an anti-quasiparticle). How to deter- 
mine M]_ is discussed in Sec. 15.11 Note that {x + y)/2 in M]_{{x + y)/2) is a macroscopic 
spacetime coordinates. It is well known that an introduction of ultraviolet renormalization 
counter terms in vacuum theory is sufficient for the out of equilibrium theory to be renor- 
malized. So we do not explicitly write down those counter terms. We assume, for simplicity, 
that the density matrix p commutes with the charge operator Q: [p, Q] = (Generalization 
to the case [p,Q] 7^ is straightforward.). 

The CTP formalism is formulated ^ by introducing an oriented closed-time path C{= 
Ci + C2) in a complex time plane. Ci runs along real xo-axis from —00 to +00 — iO^, and C2 
comes back from +00 — iO~^ to —00 — 2i0^. The real-time formalism is achieved by doubling 
every degree of freedom (f) (01,02), where 0i = 0(a;o,x) with xq G Ci is called a type-1 
field and 02 = (p{xo,x.) with xq G C2 is called a type-2 field. The classical contour action is 
written as 

J dxo J (ix£(0(x), 0^(x)) = J dxo J dyijC{x), (11) 
£^£(01,01) - £(02,05). (12) 
£ is called a hat-Lagrangian j3] . Introducing a doublet notation. 



it = (,#.;. 4), 



h=(l _°). (13) 



£0 and £ci may be written in the forms, 

£0 = -0^0;) I dS [S\x - y)dl + M\x, y)] hkv). (14) 
4i = 4{x) I A \M\x,y) - mH\x - y)] U{y). (15) 



The vertex factors may be read off from Eqs.(|7j), (fT^ and p5p . The vertex factor for a 
four-point vertex is extracted from Li^^_: 

rV = ( ^° ) = -af3. (16) 

Here, iVu = —i\ [iV22 = «A] is the vertex factor for a vertex of type-1 [type-2] fields. For a 
two-point vertex coming from £^1, we obtain 

zV;i(x, y) = I (M^{x, y) - m^6\x - y)) fg. (17) 
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In momentum space, iVd becomes 



^3, 



Throughout this paper we do not deal with initial-correlations (see, e.g., 0). 
2.2 Free quasiparticle propagator 

In this subsection, we derive the form for the free quasiparticle propagator with the help of 
quasiparticle representation [HEIIZI- 



2.2.1 Preliminary 

Following standard procedure the four kinds of propagators emerge: 

Ai2(a;,?/) = -iTr[(f)l{y)(j)i{x)p\, 
^2i{x,y) = -iTT\(f)2{x)(f)\{y)p\ , 
^22{x,y) = -iTi 



f{Mx)4iy)}p 



(19) 
(20) 
(21) 
(22) 



where p is the density matrix, and T(T) is the time-ordering (antitime-ordering) symbol, 
01, 01, 02 and 02 are the interaction-picture fields. At the end of calculation we set 0i = 02 [Ij- 
Let us introduce a matrix propagator A{x,y), where the "caret" denotes the 2x2 matrix, 
whose (ij) component is Aij{x, y). The matrix self-energy S(x, y) is defined similarly. A(a;, y) 
is defined by (c.f. Eq.dUl)) 

j <fz(5\x-z)dl+M\x,z))A{z,y) = J d^zA{x, z) (^^j\z - y) + M\z,y)^ 

= -f,5\x-y). (23) 
In the sequel, we use the short-hand notation H = F ■ G, whose "(x, y) component" is 

H{x, y) = [F- G] {x, y) = J d'zF{x, z)G{z, y). (24) 
2.2.2 Quasiparticle representation 

Motivated by equilibrium thermal field theory (ETFT), we introduce, as usual a quasi- 
particle representation for A(x,?/): 



A{x,y) 
^BAix,y) 



diag (A^, -A^) 



(25) 
(26) 
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Brix v) = ( ^^""'^^ ^^"^'^^ ^ (27) 

'^^ [Six-y) 5ix-y) + fix,y))^ 

R „\ _ ( S{x -y) + f{x,y) f{x,y) \ , 
""^^"'^^ - [ S{x-y) 6{x-y))- ^^8) 

Here A/j(Ayi) is the retarded (advanced) propagator. As will be seen in Sec M.H f{x,y) is 
related to the number density. In particular, in equilibrium thermal field theory, f{x,y) = 
f{xo - yo)5{x - y) and 

J dW^-V(x) = e(po)n.(„,)(bo|) - 0{-Po) (29) 
with n+(n_) the number density of a quasiparticle (an anti-quasiparticle). 

2.2.3 Retarded (advanced) propagator A/j(^) 
A_R(x,y) [A^(x,y)] is obtained by solving 

J d'z {5\x - z)dl + M\x, z)) AniA) {z, y) 

= I d'zAniA)ix,z) (Pj\z~y)+M'{z,y)^ 

= -6\x-y), (30) 

under the retarded (advanced) boundary condition. Wigner-transforming and taking up to 
the second-order terms in the derivative expansion, we obtain, with obvious notation, 

AniA)iX;P) ^ Ag)^)(X;P) + Ag)^)(X;P), (31) 
Ag)(X;P) ^ (Ai^^;^))^- p._^.(;,^^,^)^,(,^)o. . 



where A1^(X;po) is as in Eq. lfTm) . It should be noted that no gradient term appears. 

2.2.4 Free quasiparticle propagator A 

Wigner transformation of Eq. (f^3j) yields 

A(X;P) . ( _^ )+(A,-A,)/i, 



(33) 



+1 {Ag) + Afj} A, - 1 {{Ag) - Af, /}} i„ (34) 

i± = 



2i 
1 ±1 
±1 1 



where A^(^) = A^(^)(X; P), and / = /(X; P). {...,.. ■} and {{■ 
{{A(X;P),B(X;P)}} 



■}} are defined as 



(35) 



dX^dP'' ^P^'^X^ 



+ 



d^A 



d^B 



dX^dX^ ^P^'^P'' dP'^dP'^ dX^dX^ ' 
respectively. Straightforward computation using Eq. ljH^ yields 

{Ag)+Aj)/} = {f,P^-M\X-po)]{^f + ^f), 

'dx.dx^M\X;po)dp^dpJ - 2dlf] (Ag)' - aJ) 

+ \^8P''dx^M\X;po)dp^dx.f 

+ 2 (dx.M\X;po)) (dx^M\X;po)) dp^dpj 



(36) 



(37) 



{{^ 



(0) 



R ~ ) 



/}} 



A 



(0)^ 



A 



(0)^ 



(38) 



We may rearrange Eq.dSH) in the form A(X; P) ^ A(°) + A^^) + A^^), where k^^\i = 0, 1, 2) 
stands for the jth-order term in the derivative expansion: 



A(°)(X;P) 
A«(X;P) 
A(2)(X;P) 



Ag) 



A 



(0) 
R 



A 



(0) 





-A 



(0) 



A 



(0) 
R 



A?^) /in 



(39) 
(40) 



A 



(2) 



Ug^-A? -A?, 
-l{{A2)-A?),/}}i, 



+ (Ag)-A?)/i, 



(41) 



Eq.dSni) with Eqs.dSnil-diH) yields 
1 



2i 



2P^9x./(X; P) + dpJiX; P)dx.M\X- p,) 

-\{{Af-Afj{X-P)}} = ^. 



(42) 



Except the fact that a medium effect is taken in A^^(X;po); Eq. ()42|1 leads to a "free Boltz- 
mann equation." (see remark about after Eq. ()28|1 and Sec. 15.21 below.) 
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3 Resummation of the self-energy 



In this section, we construct and analyze the self-energy resummed propagator, and disclose 
the existence of a large contribution, which spoils the consistency of perturbative framework. 



3.1 Self-energy resummed propagator 

Let S = ?/) be a self-energy matrix. The S-resummed propagator G obeys the Dyson 
equation, 

= A + A-S-G = G-S-A + A. (43) 
Multiplying Bj^^{Bj^^) from the left (right) of both sides of Eq. pHj) . we have 

G{x, y) = ARA + ^RA-t-G = ARA + G-t-ARA, (44) 
where use has been made of Eq.fl^ and 



G{x,y) = B-^'-G-Gl' (45) 

Si? T^h 
-Sa 



S(x,y) = Bn-t-BL=( ). (46) 



Here 



T.R{x,y) = S11 + S12, S^(x,?/) = - (S22 + S12) , (47) 
Sx(a;,y) = Si2 + Si2-/-S2i-/ + S^-/-/-S^. (48) 

In obtaining Eqs. pU|) - PHj) . we have used the relation J2'^j=i^ij{x,y) = 0. Wigner trans- 
forming and taking up to the second-order terms in the derivative expansion, we obtain 

J:k{X;P) ~ zf(P) + z{/,/?eS,?}-i{{/,Si2-S2i}}, (49) 
tT^^\X;P) ^ (l + /)Si2-/S2i, (50) 

where use has been made of S^(X;P) = [T,ji(X; P)]* . T^^^ is proportional to the net pro- 
duction rate of quasiparticles (see below after Eq. ljH^ ). Eq.(jlH) with Eqs.(j2ni) t^Us 
us that G is of the form, 

with 

GR(A)ix,y) = Aij(A) + Aji(A) ■ '^R(A) ■ Gr(^a), (52) 
GKix,y) = -Gr-T.k-Ga. (53) 
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Wigner transforming and solving Eq. (j3^ . we obtain 



Gr{X;P) 



[Ga{X;P)Y C:.G^^\X-P) + G^^>{X;P), 
1 

P2 _ _^2(x; _ P) + ze(po)0+ 

gUX;P) (Ag)(X;P))"^ 
P2 - ^2(x; - S^,(X; P) + te{po)0+ 



.(2), 



GgHX;P) 



where 



-lA(°'{{E„.G2>}}-l{{A<»>E,Gi?'}} 

We may write P) in the form, 

1 



Gr{X;P) 



P2 - M^{X; po) - P) - ^(X; P) ^ 



where 



^(X;P)= [A2HX;P)Gg)(X;P)]" ^^(X;P). 
For Gx(X; P), we obtain from Eq. (jHH|) . after some manipulation, 

Gk{X;P) ~ 



l^GfGf (j^KG^i^Gf E^} + {P^ - M\ - ^a} 



+dpEK [{dx.^R + dx.M^) Gg) - (dx.^A + dx.M^) G 



—dxt^^K 



2Pn G 



'R 



dp.,T.A -2Pn G 



(54) 
(55) 

(56) 



(57) 

(58) 
(59) 



(60) 



In Eg. ()60|1 . we have kept up to the first-order term of the derivative expansion. The reason 
is as follows: In the physical-X scheme, which will be constructed in SecEl below, Sx is 
proportional to dN^ /dX^, with N~^{N~) the "physical" number density of the quasiparticle 
(anti-quasiparticle), and already is of the same order of magnitude as the first-order term of 
the derivative expansion. As has been shown in [7] and as will be briefly discussed in SecsjU 
and ini below, the bare-X scheme is equivalent to the physical-X scheme in the sense that 
they lead to the same result for the physical quantities. Then, the above statement applies 
also to the bare-X scheme. 
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We are now in a position to obtain a form for G{X;P) from Eq. ()45p . Straightforward 
but tedious manipulation yields 



where 



G{X-P) 



G^^\X; P) + P) + P) + ^^(X; P)i 



(0) 



(^(0) ^(0) ^ 



(0) 
A 



+ 



(0) 
R 



A, 



G 



(2) 



^(2) ^(2) 



+ 



(2) 



/ 



A, 



(0) 
R 



G 



(0) 

A 7 



{{gL°'-g?'./}} 



{/,p^-a<^-s4g$ 

+ {/,p2-A<2_5^^|C'( 
{{-P^ + Al^ + S«,/}}Gf 



+ 



+2 [dx.^R + dx.M^) {dx^J:R + dx^M^) dp^dpj 

+2 (dpEn - 2P^) {dp^J^R - dx.dx^f]G 
-{R-.A). 



(61) 
(62) 

(63) 



(64) 



(65) 



(66) 



Inspection of Eqs.(|UT |) -(|U ^ shows that all but GkA+ do not bring about any pathology into 
perturbation theory. We observe that Gk in Eq. ljUIH) includes the product G^I^\X; P)G^^\X; P). 
As seen from Eq. fl55|l . in the narrow width approximation, ImTiR{= —ImLA) "£(^0)0"*", 
Gr^G^^^ develops pinch singularity in a complex po-plane. Then the contribution of Gk{X; P) 
to some amplitude diverges on the energy shell, ReG~^{X] P) = (c.f. Eq. (jKK|) ). In practice, 
ImGp^{(x A^) is a small quantity, so that the contribution , although not divergent, is large. 
This invalidates the perturbative scheme. 

In the next two sections, we present two perturbative schemes, the bare-A^ scheme and 
the physical- scheme. 



4 Bare-A^ scheme 

In this section, we briefly describe the bare-A^ scheme j?]. 
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4.1 Preliminary 



In order to clarify the physical meaning of /, we compute the statistical average of current 
density, 

i 
2 



where d'^= — Taking the statistical average of Eq. ()67p . we obtain 



(67) 



(2^ 



(68) 



where Gc{X]P) = Gn + G22 = Gu + 6*21- Now we first consider the contribution from 
G'(°)(X; P) = Gf^ + Gf^, the leading part of Gc{X; P) in the derivative expansion. When 
the interaction is switched off, G^^^ reduces to A(°) = A^"^ + A22^ with M'^{X;po) = m^. The 



contribution from A^''^ with AA'^{X;pq) 



m 



(see Eqs.(JSni) with Eq.^) is 



(/(X))-(0|/(X)|0) 



(j(X))-(0|j(X)|0) 



d^p 

(27r)3 
d^P p 



[/(X;Ep,p) + {l + /(X;-Ep,p)}], 
[/(X;Ep,p)-{l + /(X;-E„p)}] 



(69) 



(27r)3Ep 

where Ep = y/p'^ + m'^, /(X; ±Ep, p) = /(X; P)\pf^=±Ep- These are to be compared with 

d^P 



(70) 



inx)) - {o\fix)\o) 



X+(X;Ep,p)-X-(X;£;p,p) 



(71) 



where = (l,p/_E'p) is the four-velocity and X"'"(X; £'p, p) [X (X;_E'p,p)] is the number 
density of the quasiparticle [anti-quasiparticle] with momentum p. Comparing Eqs. ()69|l and 
((701) with Eq.dni), we have 

/(X;po = Ep,p) = X+(X;Ep,p), /(X;po = -Ep,p) = -1 - {X; Ep, -p). (72) 

Thus, the physical number densities are obtained through computing Eq. (jHH|) . Eq. ()68|1 with 
for Gc leads to a perturbative correction to Eqs. ljU^ and (fTn)) . due to quantum and 
medium effect. This is also the case for and G^^) (c.f. Eqs.dHl - On the contrary, 

as has been discussed at the end of Sec. El the contribution from Gx(X; P) (g Gc{X; P)) is 
disastrously large. 



4.2 Bare-A^ scheme 

In this subsection, we write f^^^ for / above. The "bare" number densities N^^^^(X;pq, p) = 
eipo)f^^\X;p) and X(^)-(X; |po|,p) = -1 - e{-po)f^^\X;po,-p) obey the "free Boltz- 
mann equation", Eq. ()42|1 with f^^^ for /, which is to be solved under the given initial data 
/(^)(X°, X; P) = /(X°, X; P) with X^ the initial time. 
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The physical number densities, which are obtained from {j^{X)), Eq. ljUHj) . are functionals 
of /(i^): 

/(P^)(X;P) = 0(j9o)iV(P'^)+(X;po,p) + ^(-Po)[l-iV(P^^)-(X;bo|,-py 

= n{X;P;m. (73) 

T-C here contains large contributions mentioned above. Solving this equation for f^^\ one 
obtains 

/(B) = /(B)(X;P;[/(Ph)]). (74) 

Computation of some physical quantity yields the expression F{[f^^^]), which includes large 
contribution. Substituting the RHS of Eq. (f7^ for /*^b) one obtains the expression 

which does not include large contribution. The perturbation theory thus con- 
structed is called the "bare-iV scheme" in 



5 Physical-A^ scheme, and Gap and Boltzmann equa- 
tions 

In this section, we construct a "physical-X scheme" in terms of the number densities that 
are as close as possible to the physical number density. To this end, introducing a new 
function f{x,y), we redefine A{x,y) by Eq. pH|) with this new f{x,y) in -Bl(_r)- Then, the 
"free Boltzmann equation" ()42j) does not hold. Specification of / will be made in Sec. 15.21 
Now, in contrast to Eq.(j2SI), A is not an inverse of —f3(^6'^{x — y)dy+A4'^{x,y)^, so 

that Cq in Eq. ()14j) is not the free hat-Lagrangian. For the purpose of finding the free 
hat-Lagrangian, we first compute / d'^zf^^S^^x — z)d1 + ^A'^{x, z))K{z,y) with K{z,y) as in 
Eq.(|2i: 

d^z% {5\x - z)dl + M^[x, z)) A(z, y) 
= -6'^{x - y) 

+ r3 



O^Bl -BLd^]+M'-BL-BL-M'}-BL'-A 



{x,y), (75) 



where we have used the short-hand notation 
Eq. (ffH|) . we obtain 

j d'^z'D{x, z)A{z,y) ^ 5'^{x — y), 
V{x,y) = 



. Substituting Eq.(j2Zj) for B^ in A in 

(76) 

-f,{6\x-y)dl + M\x,y)) 
+^ J ^e-^(^-^){/,p2-A^2(X;po)}i- (77) 
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Eq. (f7njl tells us that the free hat-Lagrangian is 

4 = 1 d'y^\x)V{x,y)4>{y). (78) 
Since C'q ^ Cq, Eqs.dZH} and ((HI), the counter hat-Lagrangian Cq — C'q = tc2 appears: 

L2 = -^\x) J A / ^^e-^P<-'yh{f,P'-M'{X;po)}Aj{y). (79) 

jC.c2 yields the 2-point vertex function iVc2 = {f,P^ — -M'^{X;pQ)}A_. Thus, the theory 
contains three types of vertices, iV (Eq. (fT6|) ). iVd (Eq. lfTTj) ). and iVc2- In particular, iVci 
and iVc2 contribute to the self-energy S, 

S = S(l-P) + _ _^2(^. p^))^^ ^ ^ p2 _ _^2(^. i_. (80) 

Here S('°°p) is the contribution from the loop diagrams. 



5.1 Gap equation 



For the purpose of determining so far arbitrary mass function Ai{X;pQ), different prescrip- 
tions are available (see, e.g., iJUj.) Here we employ the pole prescription. 



Re 



0, 



po=±M±(X),p=0 



51) 



with Gr{X;P) as in Eq. ()58j) . Incorporating the contribution from the counter Lagrangian 
Cci (see Eq.dHDl)), Eq.(jHII) becomes 



M|(X) -m'- i?eSiJ°°P)(X; P) - ReJ^{X; P) 



po=±M±{X),p=0 



12) 



where J^{X;P) is as in Eq. (|59|l . Eq. (j82|l is the gap equation, from which M|(X) is deter- 
mined self consistently. It should be noted that, in Eq. (jHT|l . instead of choosing p = 0, one 
can choose p = Pc(7^ 0). 



5.2 Boltzmann equation 

We now define /, such that the number densities (c.f. Eq. lff^ ) 

iV+(X;po = ^+,p) = /(X;po = ^+,p), (83) 
iV-(X;-po = ^-,-p) = -l-/(X;po = -^-,p), (84) 
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are as close as possible to the physical- number densities. Here uj± = ujj-{X; ±p) is the energy 
of the quasiparticle (anti-quasiparticle) mode with momentum ±p: 



Re 



G-n\X;P) 

P^-rn'- i?eS^""P^(X; P) - ReJ^{X; P) 



PO=±i^±(X;±p) 
,2 D„\^(loop)/ 



0. 



Po=±'^±(^;±p) 



^5) 



As has been stressed in Sec. 14.11 Gk in Eq.ijHDI), which includes T^k (see Eq. (|H!?jl ). yields a 
large contribution to the physical number density. Then, as the determining equation for /, 
we adopt 

Ex(X;P) = 0. (86) 

It should be emphasized here that this condition is by no means unique 0. As a matter of 
fact, no large contribution emerges, as long as we impose the condition ()8fj|l on the energy 
shells and their vicinities. 

Recalling Eq. ljHUj) with Eq. (jIH|) . the condition (jHUj) turns out to 

l:t''\x■,P)-^{f,P'-m'} = 0. (87) 

Performing the derivative expansion and taking up to second-order terms, Eq. ()87|) with 
Eq. ljl^ becomes 

2P'^dx.f{X-P) - {f{X;P),M\X;po) + ReEn{X;P)} 
= f!rip)(X;P) 

+ {{/(X; P), J:^lr\X; P) - s('r)(X; P)}} . (88) 

We are now in a position to inspect the physical implication of Eq. (l88|l . Straightforward 
manipulation using Eqs. ljHHjl and (see Appendix shows that Eq. (jHH|) becomes, on the 
energy-shells, po = ±i^±, 



duj+ dN^ 



dx^ ^ ^ dXi^ dP„ 



± 



2uj± 



dXi" dp, 



dP^ ^X^' 



+ [[±N^{X- ±a;±, ±p), S_(X; P)}} 



^9) 



- PO = ±iL)± 



where v±(= ±(9u;±/(9p) is the group velocity of the ± mode with momentum ±p, Z± is the 
wave-function renormalization factor (cf. Eq. ()l()H|l ). and T is as in Eq. ljKnjl . In Eq. ljHTHl . 



1 + iV* E 



±\ v^(loop) 



±^(loop) 



"12(21) 



21(12) 



PO=±i^± 



(90) 
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E_(.Y;P) 



^(loop) 
^12 



^(loop) 
^21 



{X;P), 
<fN^ 

T2 



(91) 



dXi'dP^ dP^dX^ 



± 



d'^N^ 



dP^dP^ ^X^'^X'^ 



-n. 



± 



+ 



dPdPidX^dX^ 



-2 



dXf'dPjdP^dXi 



dXt^dX^dP^^ 



dX'dXi dXi" 



dP^dXi 



driQ 



-2 



(9X^(9P^ dX 
S^E 



SX^SX'^ dP^ 
S^E 



dP^dX'^ dP^ 



(92) 



where i and j run over 1,2,3 and Uq = dN^/dp^. r|^^ that comes from r||'^^^p^ fEq. lj^ ) is 
the net production rate of the quasiparticle (anti-quasiparticle) with momentum p(— p). In 
fact, r|^^ is the difference between the production rate and the decay rate. 
Several observations are in order. 

1) The last term on the LHS and the last term on the RHS of Eq. ()89j] are new, which 
comes from the second-order parts of the derivative expansion. 

2) The third term on the LHS of Eq.()89|) comes from the force ±duj±/dX^ due to the 
variation of the energy u)±. with X^. 

3) If the terms mentioned in 1) and 2) above are unimportant, and the renormalization of 
the wave function is neglected {Z± ~ 1), Eq.ijHHI) becomes the standard transport equation. 



9X0 



(93) 



Thus Eq. ()89j) is a generalized transport or Boltzmann equation. 

4) The last term on the LHS of Eq. ()89p is regarded as corrections to the first three terms: 



LHS of Eq. 



1 ± 



± 



dpo 
Z± 



9X0 



2uj± 
duj± 

aX^ ~ 2uj± dX^ 



+ 



2uj± 9p 



■ VxX^ 



dReJ^ 



(94) 



On the other hand the last term on the RHS of Eq. (|Hn|) is regarded as a correction to the 
production rate, which comes from the fact that the spacetime point, at which a quasi- 
(anti)particle production (decay) occurs, does not coincide with the spacetime point of 
X±(X;c2;±,±p). (cf. Eq.gHD). 

5) On the light of the above observations, one can regard Eq. ()88|1 as an off the energy-shell 
generalization of the generalized Boltzmann equation. 
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Finally, it is worth mentioning that the physical number densities computed in the 
physical- scheme contain no large terms: 

jY(ph)±(^^. ^ N^{X; P) + AN^{X; P) + . . . , (95) 
with AA^^ the perturbative correction. 



5.3 Procedure of solving the gap and Boltzmann equation 

The mass function A4'^{X;pq) and the number-density function /(X; P) are determined by 
simultaneously solving the gap (jH^ and the Boltzmann (|HH|) equations. 

As mentioned in Sec. ^ we are interested in quasiuniform systems near equilibrium and 
nonequilibrium quasistationary system. Eq. tells us that the change of along the 
(1, v-t) is of O(A^). Then, it is natural to assume that 9xm/(X; P) is also of O(A^). 

One can solve this set of equations iteratively in perturbation theory. Up to the third 
order of perturbation theory, the gap and the Boltzmann equations are dealt with to the 
gradient approximation, i.e., the last terms on the LHS of Eq. ljH^ . and on the RHS of Eq. (|55j) 
are dropped. 

1) Zeroth order of perturbation theory. 
The gap and Boltzmann equations become, 

M^^iX)=m\ (96) 
p,dfo{X;P) 

where M^^(X) and fo{X; P) are the 0(A°) = 0(1), or "bare" parts of M|(X) and /(X; P), 
respectively. Solve this Boltzmann equation under the given initial data. 

2) First order of perturbation theory. 

To 0(A), a tadpole diagram contributes to S(i°°p) = (s(i°°p)(i)) . It can readily be shown 

that (= is real and P-independent, and sfr''^^'^ = ^^21°"^'-'^ = 0- Then, 

the gap and the Boltzmann equations become 



eJ°p)(')(X)-M2^(X) 



= 0, (98) 

po=±M±{X),p=0 



and 



.r., 9fi{X;P) , g4^)(X)g(/o(X;P) + A(X;P)) ^ 

+ 9XM dP^ " ' ^ ^ 

respectively. Here Mf^{X) and /i(X; P) are the 0(A) contributions to M^(X) and /(X; P), 
respectively. Incidentally, in computing Eji.°°^'*^^''(X), we may use the leading part A(°)(X; P) 
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of the propagator A(X; P) (see, Eq-dHHl)) with M"^ Ml + Mf (= + Ml) and / ^ 
/o + /i- Eq.(pH|) is immediately solved to give 

M^^{X) = M^_{X) = (100) 

Eq.(jnni) is to be solved under the initial data /i(Xoi„, X; P) = with Xoi„ the initial time. 

3) Higher order calculation. 
At the O(A^) stage, i'^^oop)i^'i Eq.(jHHl) appears first. At the O(A^) stage, a new tadpole 
diagram appears, which is dealt with in a same manner as in the 0(A) stage. At the O(A^) 
stage, the gap and the Boltzmann equations deduced in the present paper, by taking up to 
the second-order terms of derivative expansion, enters the stage. 



6 Discussion 

Here we would like to mention a similarity between the two schemes presented here, the 
bare-iV scheme and the physical-iV scheme, and those in the ultra-violet renormalization 
schemes in quantum field theory. For simplicity, we focus on the mass renormalization. 

6.1 Structure of UV- renormalization 

"Bare" UV-renormalization scheme: 

The free Lagrangian density reads Cq = —(f)'^{x){d'^ + m%)(j){x) with rriB the bare mass. 
Computation of the physical mass mph yields mph = mpi^^mB), which includes diverging 
terms. Solving this equation for m^, we obtain uib = mBirn^Y^). Some physical quantity 
F, which we compute perturbatively, yields the expression F = F{mB), which contains, in 
general, UV-divergences. Substituting the equation tjib = mB{mpii) for ttib in F^ttib), one 
gets Fji^mph) = F{mB{mph)), which is free from UV-divergence. 

"Physical" UV-renormalization scheme: 
We introduce new free Lagrangian £[, = —(j)'^{x) (9^ + m^) 0(x) with m the renormalized 
mass. Then, the counter Lagrangian should be introduced, Cc = l^m—l^ra — 0^(3^) ["^^ ~ "^b\ 

— is determined so that the perturbatively computed physical mass is free from 
the UV-divergences. Thus, no diverging term is involved in the relation mph = mph(m). 

6.2 Structure of the two schemes presented above 

Bare-N scheme: 

The counter Lagrangian £^2 is absent and / obeys the "free Boltzmann equation." The 
physical number densities, which are the functionals of f^^\ include large contributions. 
Perturbative computation of some quantity yields the expression, which is written in terms 
of f^^^ and includes large contributions. Rewriting it in terms of the physical quantity f^'^^\ 
one obtains the large contribution free form. 
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Physical-N scheme: 

We introduce a counter Lagrangian tc2 (f7^ - which is determined so that the perturbatively 
computed physical number densities do not involve large contributions (cf. Eq. ()95|) ). 

6.3 Correspondence 

Above observation discloses the correspondence between the two schemes presented here and 
those in the UV-renormalization scheme: 
Bare scheme: 

Co{x) ^ C^ + C,i (Eq.®), 

Physical scheme: 

C',{x) ^ 4 (Eq.dZi), 
C,{x) ^ L2 (Eq-dZni)), 

m ^ /, 

mph = mph /^P") = /^P'^H/] (Eq.®), 

Absence of 
divergence ^ large contribution. 

7 Summary 

In this paper, for dealing with the out of equilibrium complex-scalar systems, we have con- 
structed two kind of perturbation theories, the bare-A^ scheme and the physical-A^ scheme, 
beyond the gradient approximation in the derivative expansion. The physical- A^ scheme has 
been constructed by imposing the condition Ex(a;, ?/) = 0, Eq. ljHTj) . This condition, which 
turns out to be a generalized Boltzmann equation on the energy-shell, serves as the deter- 
mining equation for the number density function f{X;P). At the same time, the effective 
mass M]_{X) is determined through the gap equation (jH^ . 
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A Derivation of Eq. (|HH1) 

Here we derive Eq. ljHIIjl . We define "the on-shell energy" po = ±uj±{X; ±p)(= ±uj±) through 
Eg. ()85j) ■ from which we obtain 



dRe (SS°°P^(X;P)+J^(X;P) 
dp 



-2p ± 2cu±v±Zi\ (101) 



dRe{j:T^\X;P)+nX;P)) d.^ 

where J-'{X; P) is as in Eq.(|^. v± = ±du!±/dp is the group velocity of the ± mode with 
momentum ±p, and 



1 gPe(si;°°P) + ^) 



(103) 

po=±<^± 



is the wave-function renormahzation factor. 

The determining equation for /, Eq. ()88|) . is written in the form, 



dX>^ dX^ dP, ^ dP, dX^ -^{ioop) + 8lU'^i2 ^21 n- (iU4j 
Putting Eqs. fTTrni and (|TT1^ into Eq. (fTTn|l . we obtain 



(LHSofEq.M) = 2pO^ + 2p.Vx/-^^^||^ 



ax d'P dP^ dX^' 



df dReT df dReT 



dXf' dPf, dPf, dXf' 



(105) 



Multiplying Z±/2cj± to both sides of Eq.jHH) and using Eqs. ()8Hj) and (jH3), we obtain Eq. ()89|l . 
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